DEVELOPMENT  OF  REGIONAL  TRAVEL-TIME  TABLES 
FOR  DIFFERENT  GEOTECTONIC  PROVINCES  OF  NORTHERN  EURASIA 


Victor  V.  Kirichenko  and  Yury  A.  Kraev 
Western  Services  Corporation 
Under  subcontract  to 

Science  Applications  International  Corporation 

Sponsored  by  The  Defense  Threat  Reduction  Agency 
Arms  Control  Technology  Division 
Nuclear  Treaties  Branch 

Contract  No.  DSWA01-98-C-0029 

ABSTRACT 

During  the  last  two  years,  new  regional  travel-time  tables  for  different  provinces  of  Northern  Eurasia  were 
developed  in  the  framework  of  the  Seismoacoustic  Research  for  CTBT  Monitoring  project.  This  project  is  being 
performed  by  SAIC  and  its  subcontractors:  the  Western  Services  Corporation,  the  Geophysical  Survey  of  the 
Russian  Academy  of  Sciences  (RAS)  and  Complex  Seismological  Expedition  (CSE)  of  the  Joint  Institute  of  Physics 
of  the  Earth  of  the  RAS.  One  of  the  main  objectives  of  the  project  is  to  calibrate  travel  times  for  regional  seismic 
waves  travelling  to  the  seismic  stations  of  the  Russian  Academy  of  Sciences  included  in  the  Comprehensive 
Nuclear-Test-Ban  Treaty  International  Monitoring  System  (IMS). 

The  development  of  regional  travel-time  tables  is  based  on  the  published  data  of  origin  times  and  locations  of 
nuclear  explosions  as  calibration  sources.  The  data  are  from  peaceful  nuclear  explosions  on  the  territory  of  the 
former  Soviet  Union  as  well  as  underground  nuclear  tests  at  the  Semipalatinsk  test  site,  Amchitka  Island,  and  the 
French  Test  Site  in  the  Sahara.  Values  of  arrival  times  for  major  regional  seismic  phases  (Pn/P,  Pg,  Sn/S  and  Lg) 
were  measured,  analyzed  and  collected  from  the  stations  of  the  Geophysical  Survey  and  Complex  Seismological 
Expedition. 

The  territory  of  Northern  Eurasia  was  subdivided  into  1 1  provinces  based  on  the  results  of  the  analysis  of  Pn  travel 
times  as  well  recently  published  papers  on  seismic  and  tectonic  regionalization  of  Northern  Eurasia.  At  the  22nd 
Seismic  Research  Symposium  we  present  newly  developed  travel-time  tables  for  the  following  provinces:  Central- 
East-European  territory,  Cenozoic  Folded  Regions  (Balkan  and  Carpathian  mountains,  Crimea-Caucasus-Kopet- 
Dagh,  Tyan  Shan),  Ural  folded  region,  West-Siberian  platform,  Kazakh  massif,  Altai  and  Sayan  mountains,  Siberian 
platform,  Baikal  rift  zone,  North  East  territory  and  Chukot  Peninsula.  For  all  the  developed  travel-time  tables, 
modeling  errors  were  calculated  and  their  comparison  with  the  IASPEI-9 1  tables  is  presented.  The  newly  developed 
regional  travel-time  tables  differ  substantially  from  the  IASPEI-9 1  travel-time  tables.  For  the  Pn  seismic  phase  at 
distances  of  1,500  -  2,000  km,  the  deviations  from  the  IASPEI-9 1  travel-time  tables  are  as  follows:  -8  to  -10  sec  for 
the  Central-East-European  territory,  West-Siberian  and  Siberian  platforms;  -4  to  -5  sec  for  the  Ural  folded  region, 
the  Kazakh  massif  and  the  Baikal  rift  zone  and  -2  to  +4  sec  for  the  Cenozoic  folded  regions  and  Altai  and  Sayan 
mountains.  For  the  Sn  phase  at  the  same  distances,  the  deviations  from  the  IASPEI-9 1  travel-time  tables  are, 
respectively:  -15  to  -20  sec,  -8  to- 10  sec  and  -5  to  +5  sec.  For  the  Pg  phase,  travelling  in  the  Earth's  crust,  the 
deviations  from  the  IASPEI-9 1  travel-time  tables  may  be  subdivided  into  two  groups  of  provinces  as  follows:  1) 
platforms,  Paleozoic  massifs  and  Baikal  rift  zone  and  2)  Cenozoic  folded  regions  and  Altai  and  Sayan  mountains. 
For  the  first  group  the  deviations  are  from  -3  to  -5  sec  at  the  distances  of  700  -  1,000  km.  For  the  second  group 
travel  times  are  consistent  with  the  IASPEI-9 1  travel-time  tables.  Travel  times  for  the  Lg  phase  have  large 
deviations  even  for  geotectonically  homogeneous  provinces  and  are  almost  the  same  for  the  studied  provinces.  At 
the  current  stage  of  studies  we  accepted  the  value  of  3.55  km/sec  for  the  apparent  velocity  of  Lg  phase  for  all  the 
studied  provinces.  In  general,  modeling  errors  for  the  newly  developed  travel-time  tables  are  lower  then  those  for  the 
IASPEI-9 1  tables. 

The  newly  developed  travel-time  tables  as  well  as  their  modeling  errors  were  used  for  tests  on  re-location  of 
historical  underground  nuclear  explosions  at  the  territory  of  the  former  Soviet  Union.  We  used  seismic  records  of 
the  stations  included  in  the  IMS  or  those  stations,  which  are  located  not  farther  than  300  km  from  the  IMS  stations 
within  geotectonically  homogeneous  regions.  A  comparison  between  the  mislocation  estimates  for  the  newly 
developed  travel-time  tables  and  the  IASPEI  -91  travel-time  tables  is  presented. 
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OBJECTIVE 


This  work  represents  the  results  of  a  continued  effort  to  develop  regional  travel-time  tables  for  different  geotectonic 
provinces  of  Northern  Eurasia.  The  overall  objective  is  to  calibrate  travel  times  for  regional  seismic  waves  travelling 
to  the  seismic  stations  of  the  Russian  Academy  of  Sciences  included  in  the  Comprehensive  Nuclear-Test-Ban  Treaty 
International  Monitoring  System  (IMS).  The  newly  developed  travel-time  tables  may  be  directly  used  to  improve 
seismic  event  location  as  well  as  to  test  and  to  validate  3D  travel-time  models  currently  being  developed. 


RESEARCH  ACCOMPLISHED 
Regionalization 

The  regionalization  of  the  territory  of  the  Northern  Eurasia  was  carried  out  on  the  basis  of  recent  seismicity  studies 
(Seismicity  and  Seismic  Zoning  of  Northern  Eurasia,  1993  and  1995),  geological  and  tectonic  features  of  deep 
structures  (Deep  Structure,  1991;  Pavlenko va,  1996)  as  well  as  taking  into  account  the  results  of  our  analysis  of  Pn 
travel  times.  For  presentation  of  our  results,  the  territory  of  the  Northern  Eurasia  was  subdivided  into  11  provinces: 
Central-East-European  territory,  Cenozoic  folded  regions,  Ural  folded  region,  West-Siberian  platform,  Kazakh 
massif,  Altai  and  Sayan  region,  Siberian  platform,  Baikal  rift  zone,  Amur  and  Maritime  territory,  North  East 
territory  and  Chukot  Peninsula,  Kamchatka-Kuril-Sakhalin  region.  The  polygon  vertices  for  the  specific  provinces 
are  presented  in  Table  1.  Figure  1  presents  the  regionalization  of  the  territory  of  Northern  Eurasia. 

Calibration  Events 

It  is  obvious  that  underground  nuclear  explosions  with  announced  data  on  time  origins,  locations  and  depths  are  the 
most  appropriate  calibration  events.  For  the  purposes  of  the  IMS  seismic  calibration  in  Northern  Eurasia 
underground  nuclear  explosions  conducted  at  the  territory  of  the  former  Soviet  Union  (peaceful  nuclear  explosions, 
Semipalatinsk  Test  Site  (STS)  and  Novaya  Zemlya  Test  Site),  in  the  Sahara  desert  (France)  and  on  Amchitka  Island 
(USA)  were  chosen  due  to  their  appropriate  relative  location  to  the  territory  of  the  Northern  Eurasia.  Recently 
published  data  on  peaceful  nuclear  explosions  (PNEs)  in  the  USSR  (Sultanov  et  al,  1999)  as  well  as  previously  and 
recently  published  data  on  nuclear  underground  tests  at  the  Semipalatinsk  Test  Site  (Bocharov  et  al,  1989;  Adushkin 
et  al,  1997;  Belyashova,  1999),  on  Amchitka  Island  (Springer  et  al,  1971  and  Springer  et  al,  1975)  and  in  the  Sahara 
desert  (Duclaux  et  al,  1972)  allowed  the  development  of  a  database  using  well-characterized  calibration  events. 
Calibration  events  of  the  GTO  category  were  primarily  used  for  the  development  of  regional  travel-time  tables. 

Data  on  Travel  Times 

Data  resources  were  formed  on  the  basis  of  seismic  observations  conducted  during  the  time  framework  of  1965  - 
1990  by  the  Geophysical  Survey  (GS)  of  the  RAS  and  the  Complex  Seismological  Expedition  (CSE)  of  the  Joint 
Institute  of  Physics  of  the  Earth  of  the  RAS.  Values  of  arrival  times  for  major  regional  seismic  phases  (Pn/P,  Pg, 
Sn/S  and  Lg)  were  measured,  analyzed  and  collected. 

The  total  number  of  the  GS's  stations  used  for  the  purposes  of  seismic  calibration  is  about  400.  The  total  number  of 
the  CSE's  stations  used  for  the  purposes  of  seismic  calibration  is  about  200.  Data  of  seismic  observations  of  the 
RAS's  institutions  was  added  by  the  data  of  about  300  station  in  the  Northern  Eurasia  providing  data  to  the 
International  Seismological  Center. 

Results 

The  regional  travel-time  curves  were  constructed  using  linear  regression  of  the  experimental  travel  time  data  using 
one  or  two  linear  regressions  for  the  residual  standard  deviation  minimization.  Modeling  errors  were  calculated  as 
standard  deviations  of  the  experimental  data  from  the  estimated  lines  in  a  2-degree  moving  window  with  a  50  % 
overlap.  The  basic  equation  and  the  parameters  of  the  developed  travel-time  tables  are  presented  in  the  Table  2. 
Modeling  errors  for  the  developed  travel-time  tables  are  presented  in  the  Table  3.  The  newly  developed  regional 
travel-time  tables  differ  substantially  from  the  IASPEI-91  travel-time  tables  (see  the  Abstract).  Figures  2,  3  and  4 
represents  the  comparison  of  the  Pn,  Sn,  Pg  and  Lg  reduced  travel-time  curves  for:  (1)  Central-East-European 
territory,  West-Siberian  platform,  Siberian  platform,  (2)  Ural  folded  region,  Kazakh  massif,  Baikal  rift  zone,  (3) 
Cenozoic  folded  regions,  Altai  and  Sayan  region,  North-East  territory  and  Chukot  Peninsula  and  IASPEI-91.  We 
grouped  the  aforementioned  provinces  taking  into  account  the  similarity  of  Pn  and  Sn  reduced  travel-time  tables. 


Tests  on  Relocation 

We  have  located  21  historical  underground  nuclear  explosions  at  the  STS  with  published  ground  truth  time  of 
origins  and  locations  (Bocharov  et  al,  1989)  using  the  IASPEI-91  travel-time  tables.  Then,  we  relocated  the 
aforementioned  set  of  the  explosions  using  the  newly  developed  travel-time  tables.  We  used  the  network  of  eight 
regional  seismic  stations  included  in  the  IMS  (ZAL,  MAK,  NRI,  GEYT,  BRVK,  ARU,  AAK,  TLY)  or  those 
stations,  which  are  located  not  farther  than  300  km  from  the  IMS  stations  within  geotectonically  homogeneous 
regions.  The  aforementioned  stations  are  located  at  the  range  of  6°  to  20°  from  the  underground  nuclear  explosions. 
The  defining  phases  were:  Pn,  Sn,  Pg  and  Lg.  LocSat,  provided  by  DTRA’s  Center  for  Monitoring  Research  (CMR), 
was  used  for  the  relocation  tests. 

The  testing  and  validation  metrics  are  as  follows: 

•  90%  of  events  moved  closer  to  the  GT  epicenters  with  an  average  improvement  of  8.7  km  and  a  median 
improvement  of  7.0  km; 

•  76%  of  events  moved  closer  to  the  GT  epicenters  by  20%  or  more  with  an  average  improvement  of  7.2  km  and 
median  improvement  of  7.0  km; 

•  10%  of  events  moved  away  from  GT  epicenters  with  an  average  deterioration  of  7.8  km  and  a  median 
deterioration  of  7.8  km; 

•  the  same  10%  of  events  moved  away  from  GT  epicenters  by  20%  or  more; 

•  average  mislocations  by  using  the  IASPEI-91  travel-time  tables  and  the  regional  travel-time  tables  are  18.9  km 
and  8.2  km,  respectively; 

•  median  mislocations  by  using  the  IASPEI-91  travel-time  tables  and  the  regional  travel-time  tables  are  14.9km 
and  7.9  km,  respectively; 

•  33%  of  GT  epicenters  are  within  the  calculated  90%  uncertainty  ellipses  by  using  the  IASPEI-91  travel-time 
tables  and  80%  of  GT  epicenters  are  within  the  calculated  90%  uncertainty  ellipses  by  using  the  regional  travel¬ 
time  tables; 

•  average  uncertainty  ellipse  area  is  442  km2  for  the  IASPEI-91  travel-time  tables  and  327  km2  for  the  regional 
travel-time  tables; 

•  median  uncertainty  ellipse  area  is  408  km2  for  the  IASPEI-91  travel-time  tables  and  313  km2  for  the  regional 
travel-time  tables; 

•  average  and  median  decrease  of  90%  uncertainty  ellipse  size  is  1 15  km2  and  103  km2,  respectively; 

•  average  bias  in  origin  time  estimation  is  1.1  sec  with  an  associated  standard  deviation  of  1.0  sec  (median  bias  is 
1.1  sec)  for  the  IASPEI-91  travel-time  tables  and  -0.5  sec  with  an  associated  standard  deviation  of  0.8  sec 
(median  bias  is  -0.4  sec)  for  the  regional  travel-time  tables. 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  regional  travel-time  tables  for  the  different  geotectonic  provinces  of  the  Northern  Eurasia  were  newly 
developed  and  compared  with  the  IASPEI-91  tables.  The  developed  regional  travel-time  tables  differ 
substantially  from  the  IASPEI-91  travel-time  tables. 

2.  Validation  testing  using  ground  truth  data  related  to  the  underground  nuclear  explosions  at  the  STS  shows  that 
estimates  of  event  locations  and  time  origins  are  significantly  improved  and  location  error  ellipses  are 
substantially  reduced  using  the  newly  developed  travel-time  tables. 

3.  At  the  final  stage  of  the  project  performance  we  are  planning  to  complete  the  development  of  regional  travel¬ 
time  tables  for  the  remaining  geotectonic  provinces  of  the  Northern  Eurasia,  to  clarify  the  previously  developed 
travel-time  tables  taking  into  account  still  coming  travel-time  data  from  our  subcontractors  and  to  continue 
relocation  tests  using  ground  truth  time  of  origins  and  locations  for  underground  nuclear  explosions  at  the 
territory  of  the  USSR  and  DetSSSC's  software  developed  at  the  CMR  (Xioping  Yang  et  al,  1998). 


Key  Words:  Seismic  Regional  Characterization,  Location  Calibration. 
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Table  1.  Regionalization  of  Northern  Eurasia 


Region  No. 

Region  name 

Polygon  Vertices, 

Lat.  (deg.  N)  -  Long.  (deg.  E) 

I 

Central-East-European 

72-6 

72-48 

70-48 

70-56 

territory 

48-56 

48-12 

52-12 

52-6 

II 

Cenozoic  folded  regions 

48-12 

48-62 

46-62 

46-90 

40-90 

40-78 

38-78 

38-72 

32-72 

32-42 

36-42 

36-12 

III 

Ural  folded  region 

78-48 

78-72 

74-72 

74-62 

48-62 

48-56 

70-56 

70-48 

IV 

West-Siberian  platform 

78-72 

78-102 

70-102 

70-86 

58-86 

58-82 

54-82 

54-62 

74-62 

74-72 

V 

Kazakh  massif 

54-62 

54-82 

46-82 

46-62 

VI 

Altai  and  Sayan  region 

58-82 

58-104 

46-104 

46-82 

VII 

Siberian  platform 

78-102 

78-132 

64-132 

64-138 

54-138 

54-114 

60-114 

60-104 

58-104 

58-86 

70-86 

70-102 

VIII 

Baikal  rift  zone 

60-104 

60-114 

48-114 

48-104 

IX 

Amur  and  Maritime  territory 

54-114 

54-142 

42-142 

42-114 

X 

North  East  territory 

76-132 

76-168W 

60-168W 

60-168 

And  Chukot  Peninsula 

60-156 

58-156 

58-138 

64-138 

64-132 

76-132 

XI 

Kamchatka-Kuril- 

60-156 

60-168 

52-168 

52-162 

Sakhalin  region 

48-162 

48-156 

42-156 

42-142 

54-142 

54-138 

58-138 

58-156 

100’ 


Figure  1.  Regionalization  of  Northern  Eurasia 

I.  Central-East-European  territory  II.  Cenozoic  folded  regions  III.  Ural  folded  region  IV.  West-Siberian 
platform  V.  Kazakh  massif  VI.  Altai  and  Sayan  regions  VII.  Siberian  platform  VIII.  Baikal  rift  zone 
IX.  Amur  and  Maritime  territory  X.  North  East  territory  and  Chukot  Peninsula  XI.  Kamchatka-Kuril- 
Sakhalin  region,  -  IMS  stations 


Table  2.  Parameters  of  the  Developed  Travel-Time  Curves  (Depth  =  0  km) 
Basic  Equation:  T-R/Vred  =  (A±Oa  )  -  (B  ±  oB)xR 


Phase 

Range,  km 

Vred? 

km/sec 

A 

aA 

B 

Ob 

N 

r 

SD, 

sec 

1  I.  Central-East-Euroi 

pean  territory  || 

Pn 

250-1150 

8.0 

8.15 

0.37 

0.0044 

0.0004 

75 

0.80 

0.8 

1151-2500 

8.0 

14.63 

0.39 

0.0098 

0.0002 

288 

0.93 

1.5 

Sn 

250-2500 

4.62 

15.76 

0.75 

0.0058 

0.0005 

235 

0.59 

3.4 

pg 

250-1300 

6.0 

0.64 

0.67 

0.0062 

0.0006 

75 

0.75 

1.6 

Lg 

250-2500 

3.5 

-2.25 

1.28 

0.0016 

0.0009 

167 

0.13 

5.0  i 

1  II.  Cenozoic  folded  regions  | 

Pn 

200-1800 

8.0 

8.79 

0.36 

0.0033 

0.0003 

224 

0.59 

1.9 

Sn 

250-2000 

4.62 

12.12 

1.23 

-0.0020 

0.0012 

84 

0.18 

4.4  | 

pg 

250-1200 

6.0 

1.25 

1.19 

0.0056 

0.0016 

57 

0.43 

2.4  | 

Lg 

200-2200 

3.5 

0.87 

1.63 

0.0045 

0.0016 

85 

0.29 

5.3  | 

||  III.  Ural  folded  region  || 

Pn 

250-2000 

8.0 

8.23 

0.40 

0.0048 

0.0004 

59 

0.87 

1.2 

Sn 

250-2000 

4.62 

12.86 

1.22 

-0.0002 

0.0011 

53 

0.03 

3.3 

pg 

250-1200 

6.0 

-0.60 

1.06 

0.0047 

0.0014 

35 

0.50 

1.8  | 

Lg 

200-2200 

3.5 

-2.08 

1.21 

0.0018 

0.0012 

38 

0.24 

3.1  i 

||  IV.  West-Siberian  platform  || 

Pn 

250-2300 

8.0 

8.84 

0.54 

0.0065 

0.0004 

89 

0.90 

1.7  | 

Sn 

300-2200 

4.62 

15.69 

1.27 

0.0036 

0.0009 

86 

0.38 

4.1 

pg 

250-1200 

6.0 

-0.54 

1.45 

0.0052 

0.0016 

42 

0.46 

2.3 

Lg 

250-2200 

3.5 

-1.62 

1.49 

0.0034 

0.0011 

82 

0.33 

4.8  i 

||  V.  Kazakh  massif  || 

Pn 

250-777 

8.0 

8.11 

0.39 

0.0026 

0.0006 

57 

0.51 

1.0 

778-2200 

8.0 

10.71 

0.68 

0.0061 

0.0004 

74 

0.87 

1.3 

Sn 

250-2000 

4.62 

15.96 

0.86 

0.0037 

0.0007 

82 

0.49 

2.9 

pg 

300-1200 

6.0 

0.28 

1.28 

0.0045 

0.0016 

28 

0.49 

1.9  | 

Lg 

250-2300 

3.5 

-2.24 

1.17 

0.0027 

0.0010 

86 

0.28 

4.2 

||  VI.  Altai  and  Sayan  region  || 

1  Pn 

250-1600 

8.0 

8.40 

0.62 

0.0010 

0.0008 

62 

0.23 

1.6  [ 

Sn 

300-2200 

4.62 

13.25 

1.60 

-0.0007 

0.0011 

57 

0.08 

4.2 

pg 

300-1200 

6.0 

-0.93 

3.05 

0.0026 

0.0031 

25 

0.17 

3.3  i 

Lg 

250-2200 

3.5 

-6.60 

1.47 

-0.0020 

2.1013 

99 

0.21 

4.2 

||  VII.  Siberian  platform  || 

Pn 

250-2300 

8.0 

9.83 

0.39 

0.0069 

0.0002 

206 

0.90 

1.4  | 

Sn 

600-2300 

4.62 

16.99 

1.04 

0.0067 

0.0006 

174 

0.64 

3.0 

pg 

300-1300 

6.0 

-0.90 

1.54 

0.0055 

0.0014 

29 

0.60 

1.9  | 

Lg 

250-2500 

3.5 

-3.51 

1.18 

0.00002 

0.0007 

189 

0.03 

4.0 

||  IX.  Baikal  rift  zone  || 

i  Pn 

400-2000 

8.0 

6.83 

0.59 

0.0032 

0.0006 

282 

0.87 

1.3  | 

j  Sn 

400-2000 

4.62 

12.36 

1.04 

0.0008 

0.0007 

178 

0.61 

3.3 

pg 

400-1200 

6.0 

-0.89 

1.32 

0.0037 

0.0017 

121 

0.53 

1.8  | 

Lg 

400-2500 

3.5 

-6.41 

0.59 

-0.0020 

0.0004 

245 

0.34 

4.0 

XI.  North 

East  territory  and  Chukot  Peninsula 

Pn 

350-2000 

8.0 

9.97 

0.79 

0.0044 

0.0006 

27 

0.85 

1.5 

Sn 

500-2200 

4.62 

6.62 

2.48 

-0.0050 

0.0018 

17 

0.58 

4.0 

pg 

500-1200 

6.0 

1.56 

2.31 

0.0051 

0.0023 

7 

0.69 

2.0  [ 

Lg 

400-2500 

3.5 

-0.99 

2.42 

-0.0013 

0.0014 

17 

0.23 

4.0 

Note: 


T  -  travel  time,  sec; 
N  -  data  set; 


R  -  epicentral  distance,  km;  Vred  -  reduction  velocity,  km/sec; 
r  -  correlation  coefficient;  SD  -  residual  standard  deviation,  sec. 


Table  3.  Modeling  Errors  for  the  Eurasian  Geotectonic  Provinces 


Table  3  (continued) 
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Figure  2.  Comparison  of  Pn,  Sn,  Pg  and  Lg  reduced  travel-time  curves  and  modeling  errors 

for  Central-East-European  territory  (+),  West-Siberian  platform  (*),  Siberian  platform  (  )  and 
IASPEI-91  (  ). 


TIME-R/6.0, 


Figure  3.  Comparison  of  Pn,  Sn,  Pg  and  Lg  reduced  travel-time  curves  and  modeling  errors 

for  Ural  folded  region  (+),  Kazakh  massif  (*),  Baikal  rift  zone  (  )  and  IASPEI-91  (  ). 


